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Abstract. This article gives the construction and complete classification of all 
three— dimensional spherical manifolds, and orders them by decreasing volume, in 
the context of multiconnected universe models with positive spatial curvature. It 
discusses which spherical topologies are likely to be detectable by crystallographic 
methods using three-dimensional catalogs of cosmic objects. The expected form 
^ ' of the pair separation histogram is predicted (including the location and height of 

C^") , the spikes) and is compared to computer simulations, showing that this method is 

C^") ■ stable with respect to observational uncertainties and is well suited for detecting 

C "3 ■ spherical topologies. 

o . 

PACS numbers: 98.80.-q, 04.20.-q, 02.40.Pc 

U" 1 1. Introduction 

5— 1 

(3JT), The search for the topology of our universe has focused mainly on candidate spacetimes 

with locally Euclidean or hyperbolic spatial sections. This was motivated on the one 
hand by the mathematical simplicity of three-dimensional Euclidean manifolds, and on 
the other hand by observational data which had, until recently, favored a low density 
universe. More recently, however, a combination of astrophysical and cosmological 
observations (among which the luminosity distance-redshift relation up to z ~ 1 
from type la supernovae [1], the cosmic microwave background (CMB) temperature 
anisotropies [2], gravitational lensing [3], velocity fields [4], and comoving standard 
rulers [5]) seems to indicate that the expansion of the universe is accelerating, with 
about 70% of the total energy density Oo being in the form of a dark component 
Q\ with negative pressure, usually identified with a cosmological constant term or 
a quintessence field [6]. As a consequence, the spatial sections of the universe would 
be "almost" locally flat, i.e. their curvature radius would be larger than the horizon 
radius (~ lOh^ 1 Gpc, where h is the Lemaitre-Hubble parameter in units of 100 
km/s/Mpc). Recent CMB measurements [2] report a first Doppler peak shifted by a 
few percent towards larger angular scales with respect to the peak predicted by the 
standard cold dark matter (CDM) inflationary model, thus favouring a marginally 
spherical model [7]. Indeed, under specific assumptions such as a A-CDM model and 
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a nearly scale invariant primordial power spectrum, the value of the total energy- 
density parameter is given by fl = Q mo + fl\ = 1.11+ 012 to 95% confidence [8]. 
Note that while flat models still lie well inside the 95% confidence level, the relation 
between the angular diameter distance and the acoustic peak positions in the angular 
power spectrum makes the peak positions in models with a low matter content very 
dependent on small variations of the cosmological constant. 

As a consequence of these observable facts, spherical spaceforms are of increasing 
interest for relativistic cosmology, in the framework of Friedmann-Lcmaitre solutions 
with positive spatial curvature. Due to the current constraints on the spatial curvature 
of our universe and to the rigidity theorem [9], hyperbolic topologies may be too 
large to be detectable by crystallographic methods. For example, if fl mo = 0-3 
and f^Ao = 0.6, then even in the smallest known hyperbolic topologies the distance 
from a source to its nearest topological image is more than a half of the horizon 
radius, meaning that at least one member of each pair of topological images would 
lie at a redshift too high to be easily detectable. When using statistical methods for 
detecting a hyperbolic topology, the topological signature falls to noise level as soon 
as observational uncertainties are taken into account. 

On the other hand, spherical topologies may be easily detectable, because for 
a given value of the curvature radius, spherical spaces can be as small as desired. 
There is no lower bound on their volumes, because in spherical geometry increasing a 
manifold's complexity decreases its volume, in contrast to hyperbolic geometry where 
increasing the complexity increases the volume. Thus many different spherical spaces 
would fit easily within the horizon radius, no matter how small ilo — 1 is. 

In order to clarify some misleading terminology that is currently used in the 
cosmological literature, we emphasize the distinction between spherical and closed 
universe models. A spherical universe has spatial sections with positive curvature. 
The volume of each spatial section is necessarily finite, but the spacetime can be open 
(i.e. infinite in time) if the cosmological constant is high enough. On the other hand, 
a closed universe is a model in which the scale factor reaches a finite maximum value 
before rccollapsing. It can be obtained only if the space sections are spherical and the 
cosmological constant is sufficiently low. 

As emphasized by many authors (see e.g. [10, 11, 12, 13] for reviews), the key idea 
to detect the topology in a three-dimensional data set is the topological lens effect, i.e. 
the fact that if the spatial sections of the universe have at least one characteristic size 
smaller than the spatial scale of the catalog, then different images of the same object 
shoud appear in the survey. This idea was first implemented in the crystallographic 
method [14], which uses a pair separation histogram (PSH) depicting the number of 
pairs of the catalog's objects having the same three-dimensional spatial separations in 
the universal covering space. Even if numerical simulations of this method showed the 
appearance of spikes related to characteristic distances of the fundamental polyhedron, 
we proved that sharp spikes emerge only when the holonomy group has at least one 
Clifford translation, i.e. a holonomy that translates all points the same distance [15] 
(see also [16]). Since then, various generalisations of the PSH method have been 
proposed [17, 18, 19] (see [13, 20] for a discussions of these methods) but none of them 
is fully satisfactory. The first step towards such a generalisation was to exploit the 
property that, even if there is no Clifford translation, equal distances in the universal 
covering space appear more often than just by chance. We thus reformulated the 
cosmic crystallographic method as a collecting-correlated-pairs method (CCP) [21], 
where the topological signal was enhanced by collecting all distance correlations in a 
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single index. It was proven that this signal was relevant to detect the topology. 

The goal of the present article is twofold. First we give a mathematical description 
of all spherical spaces, which have been overlooked in the literature on cosmology. This 
provides the required mathematical background for applying statistical methods to 
detect the topology: crystallographic methods using three-dimensional data sets such 
as galaxy, cluster and quasar catalogs, and methods using two-dimensional data sets 
such as temperature fluctuations in the Cosmic Microwave Background (CMB), and 
we investigate their observational signature in catalogs of discrete sources in order to 
complete our previous works [14, 15, 20, 21] on Euclidean and hyperbolic manifolds. 
We first review in § 2 the basics of cosmology and topology in spherical universes, 
including the basic relations to deal with a Fricdmann-Lcmaitre universe of constant 
curvature and the basics of cosmic topology. We then describe and classify three- 
dimensional spherical manifolds in § 3 and § 4 and explain how to use them, with 
full details gathered in Appendix A and Appendix B. One prediction of our former 
works [15, 21] was that the PSH method must exhibit spikes if there exists at least one 
Clifford translation. In § 5 we discuss the spaceforms that are likely to be detectable. 
We show that the location and height of the spikes can be predicted analytically and 
then we check our predictions numerically. We also study the effect of observational 
errors on the stability of the PSH spectra, and briefly discuss the status of the CCP 
method for spherical topologies. 

2. Cosmology in a Friedmann Lemaitre spacetime with spherical spatial 
sections 

2.1. Friedmann- Lemaitre spacetimes of constant positive curvature 

In this section we describe the cosmology and the basics of the topology of universes 
with spherical spatial sections. The local geometry of such a universe is given by a 
Fricdmann-Lcmaitre metric 

ds 2 = -dt 2 + a 2 (t) (d X 2 + sin 2 X dw 2 ) . (1) 

where a is the scale factor, t the cosmic time and dco 2 = d# 2 +sin 2 9dip 2 the infinitesimal 
solid angle, x is the (dimensionless) comoving radial distance in units of the curvature 
radius Rc of the 3-sphere S 3 . 

The 3-sphere S 3 can be embedded in four-dimensional Euclidean space by 
introducing the set of coordinates (x M ) M= o..3 related to the intrinsic coordinates 
(X, 0, tp) through (see e.g. [22]) 

x = cos x 

xi = sin x sin 9 sin tp 

X2 = sin x sin 9 cos tp 

x 3 = sin x cos 9, (2) 

with 0< x <7r, 0<#<7r and < tp < 2ir. The 3-sphere is then the submanifold of 
equation 

x^Xf, = x 2 n + x\ + xj + xj = +1, (3) 

where x^ — b lLV x v . The comoving spatial distance d between any two points x and y 
on S* 3 can be computed using the inner product x^y^. The value of this inner product 
is the same in all orthonormal coordinate systems, so without loss of generality we 
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may assume x = (1,0,0,0) and y — (cose?, sind, 0, 0), giving x^y^ = cosd. Hence, 
the comoving spatial distance between two points of comoving coordinates x and y is 
given by 

d[x, y] = arc cos [x^y^] , (4) 

The volume enclosed by a sphere of radius x is? i n units of the curvature radius, 

Vol(x) =7r(2 X -sin2 X ). (5) 

A convenient way to visualize the 3-sphere is to consider S 3 as composed of two 
solid balls in Euclidean space i? 3 , glued together along their boundaries (figure 1): 
each point of the boundary of one ball is the same as the corresponding point in the 
other ball. To represent a point of coordinates x^ in a three-dimensional space we 
consider only the coordinates (xj)j=i...3, which are located in the interior of a ball, and 
discard the nearly redundant coordinate xq. However, the two points of coordinates 
(x, 9, 4>) and (ir — x, 6, </>), corresponding respectively to the points (xq, x±, X2,xs) and 
(—xq,Xi,X2,X3) in the four-dimensional Euclidean space, have the same coordinates 
(xi, X2, X3) and we thus have to use two balls, one corresponding to < X ^ 7r/2 (i.e. 
xq > 0) and the other one to ir/2 < \ < ^ (he. xq < 0). Each ball represents half the 
space. 




Figure 1. Representation of 5 3 by two balls in R 3 glued together. Top: The 9 
and (f> coordinates are the standard ones. Bottom: The \ coordinate runs from 
at the center of one ball (the "north pole" of S 3 ) through 7r/2 at the ball's surface 
(the spherical "equator" of S 3 ) to 7r at the center of the other ball (the "south 
pole" of 5 3 ). 

All three-dimensional observations provide at least the position of an object on 
the celestial sphere and its redshift z = a/a — 1 (the value a of the scale factor today 
may be chosen arbitrarily; a natural choice is to set ao equal to the physical curvature 
radius today). To reconstruct an object's three-dimensional position (x, 8, <p) we need 
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to compute the relation between the radial coordinate \ an d the redshift z, which 
requires the law of evolution of the scale factor obtained from the Friedmann equations 

K k A. . . 

H =3 p -^ + 3 (6) 

- = -£(p + 3P) + £ (7) 

where p and P are the energy density and pressure of the cosmic fluid, A the 
cosmological constant, k = 8irG with G the Newton constant, k = +1 is the curvature 
index and H = a/a is the Hubble parameter. As a first consequence, we deduce from 
these equations that the physical curvature radius today is given by 

i?P hys = a R Co = 4r 1 = (8) 

where the density parameters are defined by 
A up 

As emphasized above, we can choose ao to be the physical curvature radius today, i.e. 
ao = RqJ s , which amounts to choosing the units on the comoving sphere such that 
Rc — 1, hence determining the value of the constant a . As long as we are dealing 
with catalogs of galaxies or clusters, we can assume that the universe is filled with a 
pressureless fluid. Integrating the radial null geodesic equation dx = At /a leads to 



X(z) = f 
Jo 



y/n mo + ^Ao - ldx 



y/n Ao + (i - n mo - n Ao )(i + x y + iyi + x y 



(10) 



2. 2. Basics of cosmic topology 

Equations (1)-(10) give the main properties that describe the local geometry, i.e. the 
geometry of the universal covering space S, independently of the topology. Indeed it 
is usually assumed that space is simply connected so that the spatial sections are the 
3-sphere S 3 . To describe the topology of these spatial sections we have to introduce 
some basic topological elements (see [10] for a review). From a topological point of 
view, it is convenient to describe a three-dimensional multi-connected manifold M by 
its fundamental polyhedron, which is convex with an even number of faces. The faces 
are identified by face-pairing isometries. The face-pairing isometries generate the 
holonomy group T, which acts without fixed points on the three-dimensional covering 
space S (see [22, 23, 24] for mathematical definitions and [10, 13] for an introduction 
to topology in the cosmological context) . The holonomy group T is isomorphic to the 
first fundamental group tti(M). 

To illustrate briefly these definitions, let us consider the particularly simple case 
of a two-dimensional flat torus T 2 . It can be constructed from a square, opposite 
edges of which are glued together. The translations taking one edge to the other are 
the face pairing isometries. The holonomy group of this space, generated by the face- 
pairing translations, is isomorphic to the group of loops on the torus tti{T 2 ) = Z x Z 
(see figure 2). 

In our case E = S 3 and its isometry group is the rotation group of the four- 
dimensional space in which it sits, i.e. G — 5*0(4). In units of the curvature radius, 
the volume of the spatial sections is given by 

Vol(S 3 /r) = ^ = |! 
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where |T| is the order of the group T, i.e. the number of elements contained in I\ 
Since the elements g of T are isometries, they satisfy 

Vx, y e E V.g e r, distfo y] = dist[ff(x), 3(2/)] (12) 

An element g of T is a Clifford translation if it translates all points the same distance, 
i.e. if 

Va;,ye£ dist[x, g(x)] = dist[y, g(y)], (13) 

The significance of these particular holonomies, which are central to the detection of 
the topology, will be explained in section 3.1. To perform computations we express 
the isometries of the holonomy group T c SO (4) as 4 x 4 matrices. To enumerate all 
spherical manifolds we will need a classification of all finite, fixed-point free subgroups 
of 50(4). This will be the purpose of sections 3 and 4. 




Figure 2. Illustration of the general topological definitions in the case of a two- 
dimensional torus (left). Its fundamental polyhedron is a square, opposite faces 
of which are identified (middle). In its universal covering space (right) the two 
translations taking the square to its nearest images generate the holonomy group. 



2.3. Spherical spaces and cosmology 

As early as 1917, de Sitter [25] distinguished the sphere 5 3 from the projective space 
P 3 (which he called elliptical space) in a cosmological context. Both spaceforms are 
finite, with (comoving) volumes 2ir 2 and it 2 , respectively. The projective space P 3 is 
constructed from the sphere 5 3 by identifying all pairs of antipodal points. The main 
difference between them is that in a sphere all straight lines starting from a given 
point reconverge at the antipodal point, whereas in projective space two straight lines 
can have at most one point in common. Thus the sphere does not satisfy Euclid's 
first axiom, while projective space does. In 5 3 the maximal distance between any two 
points is 7r, and from any given point there is only one point, the antipodal one, at 
that maximal distance. In P 3 the maximal distance is ir/2 and the set of points lying 
at maximal distance from a given point forms a two-dimensional projective plane P 2 . 
Because each pair of antipodal point points in the 3-sphere projects to a single point in 
projective space, if we adopt the 3-sphere as our space model we may be representing 
the physical world in duplicate. For that reason, de Sitter claimed that P 3 was really 
the simplest case, and that it was preferable to adopt this in a cosmological context. 
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Einstein did not share this opinion, and argued that the simple connectivity of S 3 
was physically preferable [26]. Eddington [27], Friedmann [28] and Lemaitre [29] also 
referred to projective space as a more physical alternative to S 3 . 

Narlikar and Seshadri [30] examined the conditions under which ghost images of 
celestial objects may be visible in a Friedmann-Lemaitre model with the topology 
of projective space. De Sitter had correctly remarked that the most remote points 
probably lie beyond the horizon, so that the antipodal point, if any, would remain 
unobservable. Indeed, given the present data, it will be impossible to distinguish S 3 
from P 3 observationally if Q — 1 <C 1, because our horizon radius is too small. 

None of these authors mentioned other multi-connected spaces in a cosmological 
context. This was first done by Ellis [31] and Gott [32]. More recently, lens spaces 
have been investigated in the framework of quantum gravity models [33] and in terms 
of their detectability [34]. Nevertheless, the literature on multi-connected spherical 
cosmologies is underdeveloped; here we aim to fill the gap. 



3. The Mathematics of Spherical Spaces I: Classification of S 3 subgroups 

This section describes the classification of all spherical 3-manifolds and develops an 
intuitive understanding of their topology and geometry. Of particular observational 
relevance, we will see which spherical 3-manifolds have Clifford translations in their 
holonomy groups and which do not. Thrclfall and Seifert [35] gave the first complete 
classification of these spherical 3-manifolds. Our approach borrows heavily from 
theirs, while also making use of quaternions as in [36]. 

As emphasized in the previous section, technically what we will need in order 
to perform any computation is the form of the holonomy transformations g as 4 x 4 
matrices in SO(A), and an enumeration of all finite subgroups T C 50(4). We will 
enumerate the finite subgroups of 5*0(4) (see sections 3.1 and 4) in terms of the simpler 
enumeration of finite subgroups of SO(3) (see section 3.2). The connection between 
S*0(4) and SO(3) will use quaternions (see Appendix A and Appendix B). 



3.1. Generalities 



Our starting point is the fact that for each isometry g G 0(n + 1) there is a basis of 
R n+1 relative to which the matrix of g has the form 

/ +1 ••• \ 

o '•• 

+1 



(14) 



cos a 
sin a 



- sin a 
cos a 



V 



If g is a holonomy transformation of an n-manifold, then g has no fixed points and its 
matrix has no +1 terms on the diagonal. Furthermore, each pair of —1 terms may be 
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rewritten as a sine-cosine block with a — ir. Thus when n = 3 the matrix takes the 
form 

\ 



M(9, 



/ cos 9 


— sm9 








sin 9 


cos 9 














COS (j) 


— sin 


V o 





sin <j) 


cos<? 



(15) 



An immediate consequence of this decomposition is that every spherical 3-manifold 
is orientable. Indeed all odd-dimensional spherical manifolds must be orientable for 
this same reason. In even dimensions the only spherical manifolds are the n-sphere 
S n (which is orientable) and the n-dimensional projective space P n (which is non 
orientable) . 

If we replace the static matrix M{9, (f>) with the time-dependent matrix M{9t, (pt), 
we generate a flow on S 3 , i.e. to each point x e S 3 we associate the flow line 
x(t) = M(9t, 4>t)x. This flow is most beautiful in the special case 9 = ±</>. In this 
special case all flow lines are geodesies (great circles), and the flow is homogeneous in 
the sense that there is an isometry of S 3 taking any flow line to any other flow line. 
The matrix M(9,±9) defines a Clifford translation because it translates all points 
the same distance (see equation 13). We further distinguish two families of Clifford 
translations. When 9 = <p the flow lines spiral clockwise around one another, and the 
Clifford translation is considered right-handed whereas when 9 — —<j) the flow lines 
spiral anticlockwise around one another, and the Clifford translation is considered 
left-handed. 

Every isometry M(6, (f>) € 50(4) is the product of a right-handed Clifford 
translation M(a, a) and a left-handed Clifford translation M((3, —(3) as 

M(9, 4>) = M(a, a) M(j3, -(3) = M(f3, -/?) M(a, a) (16) 

where a = (9 + <fi)/2 and (3 = (9 — (j>)/2 and the order of the factors makes no 
difference. This factorization is unique up to simultaneously multiplying both factors 
by -1. Moreover every right-handed Clifford translation commutes with every left- 
handed one, because there is always a coordinate system that simultaneously brings 
both into their canonical form (15). 

Just as the unit circle S 1 enjoys a group structure as the set S 1 of complex 
numbers of unit length, the 3 sphere S 3 enjoys a group structure as the set <S 3 of 
quaternions of unit length (see Appendix A for details). Each right-handed Clifford 
translation corresponds to left multiplication by a unit length quaternion (q — > xq), 
so the group of all right-handed Clifford translations is isomorphic to the group S 3 of 
unit length quaternions, and similarly for the left-handed Clifford translations, which 
correspond to right multiplication (q — > qx). It follows that SO (A) is isomorphic 
to <S 3 x iS 3 /{±(l, 1)}, where 1 is the identity quaternion, so the classification of all 
subgroups of SO (4) can be deduced from the classification of all subgroups of S 3 . 

Classifying all finite subgroups of S 3 seems difficult at first, but luckily it reduces 
to a simpler problem. The key is to consider the action of the quaternions by 
conjugation. That is, for each unit length quaternion x <E <S 3 , consider the isometry 
p x that sends each quaternion q to xqx -1 

f S 3 — > S 3 
P * : { q~Mq)=xqx-i ' 

The isometry p x fixes the identity quaternion 1, so in effect its action is confined to the 
equatorial 2-sphere spanned by the remaining basis quaternions (i, j, k) [see Appendix 
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A for details and definitions concerning quaternions]. Thus, by restricting our 
attention to the equatorial 2-sphere, we get an isometry 



In other words, each x £ S 3 defines an element p x G 50(3), and the mapping 



is a homomorphism from S 3 to 5*0(3). To classify all subgroups of S 3 we must first 
know the finite subgroups of 5*0(3). 

3.2. Finite Subgroups of 50(3) 

The finite subgroups of 50(3) are just the finite rotation groups of a 2-sphere, which 
are known to be precisely the following: 

• The cyclic groups Z n of order n, generated by a rotation through an angle 27r/n 
about some axis. 

• The dihedral groups D m of order 2m, generated by a rotation through an angle 
2n/m about some axis as well as a half turn about some perpendicular axis. 

• The tetrahedral group T of order 12 consisting of all orientation-preserving 
symmetries of a regular tetrahedron. 

• The octahedral group O of order 24 consisting of all orientation-preserving 
symmetries of a regular octahedron. 

• The icosahedral group I of order 60 consisting of all orientation-preserving 
symmetries of a regular icosahedron. 

If the homomorphism p : S 3 —* SO (3) were an isomorphism, the above list 
would give the finite subgroups of S 3 directly. We are not quite that lucky, but 
almost: the homomorphism p is two-to-one. It is easy to see that p x = p_ x because 
xqx 1 = (— x)q(— x) _1 for all q. There are no other redundancies, so the kernel of p 
is 



Let T be a finite subgroup of S 3 and consider separately the cases that T does or 
does not contain —1. 

• Case 1: If T does not contain — 1, then p maps T one-to-one onto its image in 
50(3), and T is isomorphic to one of the groups in the above list. Moreover, 
because we have excluded —1, and S 3 contains no other elements of order 2, we 
know that T contains no elements of order 2. The only groups on the above list 
without elements of order 2 are the cyclic groups Z n of odd order. Thus only 
cyclic groups of odd order may map isomorphically from S 3 into 50(3), and it is 
easy to check that they all do, for example by choosing the generator of Z n to be 
the quaternion cos(27r/n) 1 + sin(27r/n) i. 

• Case 2: If T contains — 1, then p maps T two-to-one onto its image in 50(3), 
and T is a two-fold cover of one of the groups in the above list. Conversely, every 
group on the list lifts to a group T C S 3 . The construction of the two-fold cover is 
trivially easy: just take the preimage of the group under the action of p. The result 
is called the binary cyclic, binary dihedral, binary tetrahedral, binary octahedral, 



Px : 5 2 5 2 



(18) 




(19) 



Ker(p) = {±1}. 



(20) 
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or binary icosahedral group. A "binary cyclic group" is just a cyclic group of twice 
the order, so all even order cyclic groups occur in this fashion. The remaining 
binary groups are not merely the product of the original polyhedral group with 
a Z 2 factor, nor are they isomorphic to the so-called extended groups which 
include the orientation-reversing as well as the orientation-preserving symmetries 
of the given polyhedron, but are something completely newf. Note that the 
plain dihedral, tetrahedral, octahedral, and icosahedral groups do not occur as 
subgroups as of S 3 - only their binary covers do. 

3.3. Finite Subgroups of S 3 

Combining the results of the two previous cases, we get the complete classification of 
finite subgroups of the group S 3 of unit length quaternions as follows: 

• The cyclic groups Z n of order n. 

• The binary dihedral groups D* m of order 4m, m > 2. 

• The binary tetrahedral group T* of order 24. 

• The binary octahedral group O* of order 48. 

• The binary icosahedral group I* of order 120. 

4. The Mathematics of Spherical Spaces II: Classification of spherical 
spaceforms 

There are three categories of spherical 3-manifolds. The single action manifolds are 
those for which a subgroup R of S 3 acts as pure right-handed Clifford translations. 
The double action manifolds arc those for which subgroups R and L of S 3 act 
simultaneously as right- and left-handed Clifford translations, and every element of R 
occurs with every element of L. The linked action manifolds are similar to the double 
action manifolds, except that each element of R occurs with only some of the elements 
of L. 

After introducing some definitions, we give the classifications of single action 
manifolds (§ 4.1), double action manifolds (§ 4.2) and linked action manifolds (§ 4.3) 
and present in § 4.4 a summary of these classifications. 

We define a lens space L(p, q) by identifying the lower surface of a lens-shaped 
solid to the upper surface with a 2-xq/p rotation (see figure 3), for relatively prime 
integers p and q with < q < p. Furthermore, we may restrict our attention to 
< q < p/2 because for values of q in the range p/2 < q < p the twist 2nq/p is the 
same as —2ir(p—q)/p, thus L(p, q) is the mirror image of L(p,p — q). When the lens is 
drawn in Euclidean space its faces are convex, but when it is realized in the 3-sphcrc 
its faces lie on great 2-spheres, filling a hemisphere of each. Exactly p copies of the 
lens tile the universal cover S 3 , just as the 2-dimensional surface of an orange may 
be tiled with p sections of orange peel, each of which is a bigon with straight sides 
meeting at the poles. Two lens spaces L(p, q) and L(p' , q') are homeomorphic if and 
only if p = p' and either q = ±q'(modp) or qq' = ±l(modp). 

A cyclic group Z n may have several different realizations as holonomy groups 
r C 50(4). For example, the lens spaces L(5, 1) and L(5,2) are nonhomcomorphic 
manifolds, even though their holonomy groups are both isomorphic to Z 5 . For 

X The binary dihedral group is isomorphic to the plain cyclic group Z4. 



Topological Lensing in Spherical Spaces 11 

noncyclic groups, the realization as a holonomy group T C SO (4) is unique up to 
an orthonormal change of basis, and thus the resulting manifold is unique. 




Figure 3. Construction of a lens space L(p,q). 

4.1. Single Action Spherical 3-Manifolds 

The finite subgroups of S 3 give the single action manifolds directly, which are thus 
the simplest class of spherical 3-manifolds. They are all given as follows: 

• Each cyclic group Z n gives a lens space L(n, 1), whose fundamental domain is a 
lens shaped solid, n of which tile the 3-sphere. 

• Each binary dihedral group gives a prism manifold, whose fundamental 
domain is a 2m-sided prism, Am of which tile the 3-sphere. 

• The binary tetrahedral group T* gives the octahedral space, whose fundamental 
domain is a regular octahedron, 24 of which tile the 3-sphere in the pattern of a 
regular 24-cell. 

• The binary octahedral group O* gives the truncated cube space, whose funda- 
mental domain is a truncated cube, 48 of which tile the 3-sphere. 

• The binary icosahedral group I* gives the Poincare dodecahedral space, whose 
fundamental domain is a regular dodecahedron, 120 of which tile the 3-sphere 
in the pattern of a regular 120-cell. Poincare discovered this manifold in a 
purely topological context, as the first example of a multiply connected homology 
sphere [37]. A quarter century later Weber and Seifert glued opposite faces of 
a dodecahedron and showed that the resulting manifold was homeomorphic to 
Poincare's homology sphere [38]. 

In figure 4, we depict the fundamental domains for the binary tetrahedral group 
T* , binary octahedral group O* and binary icosahedral group I*. The fundamental 
polyhedron for the lens space L(n, 1) can be constructed by following the example 
presented in figure 3. Finally, the fundamental domain of the prism manifold generated 
by the binary dihedral group D§ is shown in figure 5. 

To finish, let us emphasize that all single action manifolds are globally homoge- 
neous, in the sense that there is an isometry h of the manifold taking any point x to 
any other point y. If the manifold's holonomy group is realized as left multiplication by 
a group r = {<7i} of quaternions, then the isometry h is realized as right multiplication 
by x~ lr y. To check that h is well-defined on the quotient manifold S 3 /V, note that h 
takes any point giX equivalent to i to a point gix(x~ 1 y) = giy equivalent to y, thus 
respecting equivalence classes of points. 

4-2. Double action spherical 3-manifolds 

The double action spherical 3-manifolds are obtained by letting one finite subgroup 
R C S 3 act as right-handed Clifford translations (equivalent to left multiplication 
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Figure 4. Fundamental domains for three single action 3-manifolds. From left to 
right, the regular octahedron, the truncated cube and the regular dodecahedron 
which respectively correspond to the spaces generated by the binary tetrahedral 
group T* , the binary octahedral group O* and the binary icosahedral group I*. 




Figure 5. An example of fundamental domain for a prism manifold. This 10— 
sided prism is the fundamental polyhedron of the space generated by the binary 
dihedral group Dg, of order 20. 

of quaternions) while a different finite subgroup L C S 3 simultaneously acts as left- 
handed Clifford translations (equivalent to right multiplication of quaternions). A 
priori any two subgroups of <S 3 could be used. However, if an element M(9,9) of R 
and an element M(<fi,—(j>) of L share the same translation distance \9\ = \<f>\, then 
their composition, which equals M(26, 0) or M(0, 20), will have fixed points unless 
9 = <f> = ir. Allowing fixed points would take us into the realm of orbifolds, which 
most cosmologists consider unphysical to describe our universe, so we do not consider 
them here. In practice this means that the groups R and L cannot contain elements of 
the same order, with the possible exception of ±1. The binary dihedral, tetrahedral, 
octahedral, and icosahedral groups all contain elements of order 4, and so cannot be 
paired with one another. 

Thus either R or L must be cyclic. Without loss of generality we may assume 
that L is cyclic. If R is also cyclic, then L and R cannot both contain elements of 
order 4, and so we may assume that it is L that has no order 4 elements. If R is 
not cyclic, then R is a binary polyhedral group, and again L can have no elements of 
order 4. Thus either L — Z n or L = Z 2n , with n odd. Recalling that the group Z n 
consists of the powers {q l }o<i<n of a quaternion q of order n, it's convenient to think 
of the group Zi n as {q*}o<i<n U {— q l }a<i< n - If R contains —1, then nothing is gained 
by including the {— q 1 } in L, because each possible element (r)(— I) already occurs 
as (— r)(l). In other words, L = Z n and L = Z 2n produce the same result, the only 
difference being that with L = Zi n each element of the resulting group is generated 
twice, once as (r)(l) and once as (— r)(— I). If R docs not contain —1, then it is cyclic 
of odd order and we swap the roles of R and L. Either way, we may assume L is cyclic 
of odd order. The double action spherical 3-manifolds are therefore the following: 

• R — Z m and L — Z n , with m and n relatively prime, always yields a lens space 
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L(mn,q). However, not all lens spaces arise in this way. 

• R = and L = Z n , with gcd(4m, n) — 1, yields an n-fold quotient of a prism 
manifold that is simultaneously a 4m-fold quotient of the lens space L(n, 1). 

• R = T* and L = Z n> with gcd(24,n) = 1, yields an n-fold quotient of the 
octahedral space that is simultaneously a 24-fold quotient of the lens space 
L{n,l). 

• R = O* and L — Z n , with gcd(48, n) — 1, yields an n-fold quotient of the 
truncated cube space that is simultaneously a 48-fold quotient of the lens space 
L{n,l). 

• R = I* and L = Z„, with gcd(120,n) — 1, yields an n-fold quotient of the 
Poincare dodecahedral space that is simultaneously a 120-fold quotient of the 
lens space L(n, 1). 

In hgure 6, we present the fundamental domain of the double action manifold 
generated by the binary octahedral group R — O* and the cyclic group L — of 
orders 48 and 5 respectively. 




Figure 6. A fundamental domain for the double action manifold of order 240 
generated by the binary octahedral group R = O* and the cyclic group L = Z$. 



4-3. Linked Action Spherical 3-Manifolds 

The third and final way to construct spherical 3-manifolds is to choose groups R and 
L as before, but allow each element r € R to pair with a restricted class of elements 
I G L, being careful to exclude combinations of r and I that would create fixed points. 
The following linked action manifolds arise in this way. 

• R and L are the cyclic groups respectively generated by 

r = M ( P + q + 1 2,, P + q+ \n 
\ 2p ' 2p 



and 

i=m[P- 

2p ' 2p 

with < q < p and gcd(p, q) = 1. The generator r is linked to the generator I, and 
their powers are linked accordingly. This yields the lens space L(p,q), generated 
by rl = M (2ir/p, 2nq/p). Note that when p + q is even, each element of L(p, q) 
is produced twice, once as r k l k and once as r P+ k lP+ k — r p r k l p l k = (—r k )(—l k ). 

R = T* and L = Zg n with n odd. The plain (not binary) tetrahedral group T 
contains a normal subgroup H ~ Di consisting of the three half turns plus the 
identity. Each element r in the binary tetrahedral group T* is assigned an index 
0, 1, or 2 according to the coset of H in which its projection p r e 5*0(3) lies. 
Each element I in Zg n is assigned an index 0, 1, or 2 equal to its residue modulo 3. 
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An element r is linked to an element I il and only if their indices are equal. This 
yields a holonomy group with only a third as many elements as the full double 
action group would have, and avoids elements with fixed points. 
• R = D* m and L = Zg, n with gcd(m, 8n) = 1. Each element r in D* m is assigned 
an index or 1 according to whether its projection p r G SO(3) lies in the cyclic 
part of the plain D m or not. Each element I in Z$ n is assigned an index or 1 
equal to its residue modulo 2. An element r is linked to an element I if and only 
if their indices are equal. This yields a holonomy group with only half as many 
elements as the full double action group would have, and avoids elements with 
fixed points. Note that because R and L both contain —1, each element in the 
group is produced twice, once as rl and once as (— r)(— I). 

As an example, figure 7 shows the fundamental polyhedron of the linked action 
manifold generated by the binary tetrahedral group R = T* and the cyclic group 
L = Zg, of orders 24 and 9 respectively. 




Figure 7. Fundamental domain of the linked action manifold of order 72 
generated by the binary tetrahedral group R = T* and the cyclic group L = Z9. 

4-4- Summary of classification 

The preceding sections constructed all three-dimensional spherical manifolds and clas- 
sified them in three families. To sum up, figure 8 organizes these manifolds by de- 
creasing volume. Since none homeomorphic lens spaces can have isomorphic holonomy 
group, figure 9 charts this special case of lens spaces. 

On Figure 10, we show the number of distinct spherical 3-manifolds of a given 
order. The vast majority of these manifolds are lens spaces. This does not necessarily 
mean that a spherical universe is "more likely" to be a lens space. It does, however, 
reflect the fact that there is a free parameter governing the amount of twist in the 
construction of a lens space of order n. The remaining groups of order n have no 
free parameters in their construction. The difference arises because each lens space 
is generated by a single element and is therefore subject to minimal consistency 
constraints; the remaining manifolds, which each requires at least two generators, 
are subject to stronger consistency conditions. 

In the construction of a lens space of order n, there are roughly n choices for 
the amount of twist (the exact number of choices varies with n, because the twist 
parameter must be relatively prime to n), so the number of distinct spherical manifolds 
of order exactly n grows linearly with n (see Figure 10, left plot). Thus the total 
number of spherical 3-manifolds of order n or less is the sum of an arithmetic series, 
and therefore grows quadratically with n (see Figure 10, right plot). 
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Z 8 x Z 27 D* x Z 27 
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■^240 Dq 


Z 8 o X Z3 Z4S X Z5 Zi6 X Z15 
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O* x Z 5 


D20 x Z3 Z240 
Dl x Z 48 D* 15 x Z 16 



Figure 8. Classification of groups generating single, double and linked action 
spherical 3-manifolds. The first column gives the order |T| of the holonomy group, 
i.e. along each row the volume of space is 2tt 2 / 1 T| . The simply connected 3-sphcrc 
is S 3 /Zi, and the projective space P 3 is S 3 /Z2- Since some cyclic groups Z n have 
different realizations as single action, double action, and linked action manifolds 
they may appear more than once on a given line. For example, the lens spaces 
L(5, 1) and 1/(5, 2) are nonhomcomorphic manifolds, even though their holonomy 
groups arc both isomorphic to Z5. The double action groups Z m X Z n ~ Z mn 
all yield lens spaces, as do the single action groups Z mn . For instance 1/(12, 1) 
is the single action manifold generated by Z12 while L(12, 5) is the double action 
manifold generated by Z?, X Z4 ~ Zi2- For noncyclic groups, the realization as a 
holonomy group is unique, and thus the resulting manifold is unique, named after 
its holonomy group. 

5. Crystallographic Simulations 

As we emphasized, the PSH method applies as long as the holonomy group has at 
least one Clifford translation. Thus we use the PSH method for most of the spherical 
manifolds but we will need the CCP method in certain exceptional cases. 

Section 5.1 determines the radius Xmax of the observable portion of the covering 3— 
sphere in units of the curvature radius, as a function of the cosmological parameters ft m 
and Oa and a redshift cutoff z c . For plausible values of x max , Section 5.2 determines 
which topologies are likely to be detectable. Section 5.3 explains the expected form of 
the Pair Separation Histogram first in the 3-sphere, and then predicts the location and 
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Figure 9. Classification of lens spaces. In this chart each lens space occurs only 
once in the first valid column (e.g. if the lens space is a single action manifold 
we ignore the trivial expression of it as a double or linked action manifold). This 
also takes into account the various equivalences so, for example L(7, 2) appears 
while L(7, 3) does not because L(7, 3) = L(7, 2). 

height of the spikes in a multiply connected spherical universe. Computer simulations 
(§ 5.4) confirm these predictions. In § 5.5, we briefly recall the applicability of the 
CCP method. 

5.1. Observational prospects 

Observations indicate that |Q TOo + Cl\ — 1 1 is at most about 0.1, and this value fixes 
the physical curvature scale i?J^ hys from (8). To detect topology one requires that 
the size of the manifold be smaller than the diameter of the observable universe 
in at least one direction. Equation (10) gives the distance \ m radians from the 
observer to a source at redshift z. Figure 11 uses equation (10) to plot the maximal 
radial distance Xmax that is accessible in a catalog of sources extending to redshift 
z c , as a function of the cosmological parameters VL mo and Oa . This maximal radial 
distance Xmax may also be used to compute the volume of the observable universe 
vol(x max ) = 7r(2x m ax - sin2x m ax) ^ f^Xmax in curvature radius units and compare 
it to the total volume 2-7T 2 of the 3-sphere. 

In practice one requires that shortest distance between topological images be less 
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Figure 10. Left: The number of distinct spherical 3— manifolds of order exactly 
|r|. The vast majority of these manifolds are lens spaces for which there are 
roughly |T| choices for the amount of twist in their construction. So, the number 
of distinct spherical manifolds of order exactly |T| grows linearly with Right: 
The total number of spherical 3— manifolds of order |T| or less is the sum of an 
arithmetic series, and therefore grows quadratically with 



than the effective cutoff radius Xmax. For example, a cyclic group Z n will satisfy this 
criterion if and only if 

< Xmax- (21) 

n 

Figure 12 shows the implications of this equation. In principle one could detect 
topology even if the shortest translations were almost as large as the diameter of 
the observable universe, i.e. 2x m ax, but when using statistical methods the signal 
would become too weak. 



5.2. Geometrical expectations 

In this section we determine which topologies are likely to be detectable for plausible 
values of Xmax (see also [34] for detectability of lens spaces) . 

In a single action manifold every holonomy is a Clifford transformation, so in 
principle the PSH detects the entire holonomy group. In practice we see only a small 
portion of the covering 3-sphere. For example, with f2 mo = 0.35, JIao = 0.75, and 
a redshift cutoff of z = 3, we see a ball of radius Xmax — 1/2 (see figure 11). With 
this horizon size the binary tetrahedral, binary octahedral, and binary icosahedral 
groups would be hard to detect because their shortest translations distances are 2tt/6, 
27r/8, and 27r/10, respectively. If, however, we extend the redshift cutoff to z = 1000, 
then the horizon radius expands to \ ~ 1, putting the binary tetrahedral, binary 
octahedral, and binary icosahedral groups within the range of CMB methods. 

The cyclic groups Z n and the binary dihedral groups D* m are much more amenable 
to detection, because the shortest translation distances, 2ir/n and 27r/2m respectively, 
can be arbitrarily small for sufficiently large n and m. For the example given above 
with Xmax — 1/2, a cyclic group Z n will be detectable if n > 27r/(l/2) ~ 12, and 
similarly a binary dihedral group will be detectable if 2m > 12. On the other hand, 



Topological Lensing in Spherical Spaces 



18 




Figure 11. The maximum radial distance Xmax accessible in a catalog of depth 
z c = 5 as a function of the cosmological parameters. 

the lack of obvious periodicity on a scale of 1/20 the horizon radius§ implies that 
n and 2m may not exceed 27r/(l/20) ~ 120. In conclusion, the orders of the cyclic 
groups Z n and the binary dihedral groups that are likely to be detectable must 
lie in the ranges 

12 ~ n ~ 120 and 12 - 2m - 120. 

The holonomy group of a double action manifold (recall § 4.2) consists of the products 
rl, as r ranges over a group R of right-handed Clifford translations and I ranges over 
a group L of left-handed Clifford translations. Because only pure Clifford translations 
generate spikes, the PSH detects a product rl if and only ifr = ±lorZ = ±l. Thus 
the set of spikes in the PSH for the double action manifold defined by R and L will be 
the union of the spikes for the two single action manifolds defined by R and L alone, 
or possibly R x {±1} or L x {±1}. For example, if R — D\ and L = Z 3 , then we 
get the union of the spikes for with the spikes for Z$ x {±1} = Z§ (sec figure 14). 
The conditions for observability thus reduce to the conditions described for a single 
action manifold. That is, a double action manifold is most easily detectable if one of 
its factors (possibly extended by —1) is a cyclic group Z n or a binary dihedral group 

§ Note that this upper bound is very approximate and that a detailed study will be required. 
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Figure 12. The redshift of the first topological image as a function of the order 
of the cyclic group Z n when Qq = 0.35 and £7a = 0.75. 



D* m , with 

12 ~ n ~ 120 and 12 - 2m - 120. 

Happily the shortest translations are all Clifford translations. 

In a linked action manifold (see § 4.3) not every element of R occurs with every 
element of L. Let R' be the subgroup consisting of those elements r E R that occur 
with ±1 of L, and let V be the subgroup consisting of those elements I £ L that occur 
with ±1 of R. With this notation, the Clifford translations are the union of R' and 
L', possibly extended by —1. 

Linked action manifolds are the most difficult to detect observationally, because 
the groups R' and L' may be much smaller than R and L. Indeed for some linked 
action manifolds the groups R 1 and L' are trivial. For example, the lens space L(5, 2) 
is obtained as follows. Let R = Z$ and L = Z§, and let the preferred generator r of 
R be a right-handed Clifford translation through an angle 2n/5 while the preferred 
generator I of L is a left-handed Clifford translation through an angle 4w/5. Link the 
generator r to the generator I, and link their powers accordingly. The identity element 
of R (resp. L) occurs only with the identity element of L (resp. R), so both R' and L' 
are trivial. In this case the PSH contains no spikes and we must use the CCP method 
or CMB methods instead. 

We emphasize that the detectability of a given topology depends not on the or- 
der of the group, but on the order of the cyclic factor (s). Given that we can see to 
a distance of about 27r/12 in comoving coordinates (see discussion above), the lens 
space L(17, 1), a single action manifold of order 17, would be detectable. However, 
the lens space L(99, 10), a double action manifold of order 99 generated by Zg x Zn, 
would probably not be detectable, because its shortest Clifford translations have dis- 
tance 27r/ll. The only way we might detect L(99, 10) would be if the Milky Way, 
just by chance, happened to lie on a screw axis, where the flow lines of the right- 
and left-handed factors coincide. On the other hand, a double action manifold like 
S 3 /(Zgr x Zgs) probably is detectable. It has Clifford translations of order 97 and 
98, so crystallographic methods should in principle work, and it does not violate the 
observed lack of periodicity at 1/20 the horizon scale because its screw axis (along 
which the minimum translation distance is 2n/(97 * 98) = 27r/9506) probably comes 
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nowhere near the Milky Way. 

Putting the results all together and using the complete classification of spherical 
spaces, we estimate that a few thousands of potentially observable topologies are 
consistent with current data. This number is reasonably low in view of search strategies 
for the detection of the shape of space through crystallogaphic or CMB methods. We 
note that similar lower bounds were recently obtained in [34] in the particular case of 
lens spaces, and this analysis does not give upper bounds and is less general than the 
one presented here. 

5.3. PSH spectra 

Let us now discuss the form of ideal PSH spectra. We start by considering the expected 
PSH spectrum in a simply-connected spherical space. As shown in [39], it is obtained 
by computing the probablility V(Xa,Xl)^Xl that two points in a ball of radius Xa are 
separated by a distance between xi and Xi + ^Xi 



valid for all Xa <= [0, 7r] and xi <= [0, min(2x a , ir)] and in which 9 stands for the 
Heavyside function. 

As seen in the sample PSH spectra of figure 14, the PSH of the simply-connected 
space S 3 provides the background contribution of the Poisson distribution over which 
the spikes of the multi-connected space will appear. 

Each Clifford translation g of the holonomy group will generate a spike in the 
PSH spectrum. A Clifford translation takes the form M(6, 9) or M(9, -9) (sec § 3.1) 
and so its corresonding spike occurs at x = The number of group elements sharing 
the same value of \0\, with 9 normalized to the inverval [— n, ir], is the multiplicity 
mult(#) of the translation distance 9. 

To compute the amplitude of the spike, we first consider the case that there is 
only one source in the fundamental polyhedron. In that case, the number of images in 
the covering 3-sphere equals the order of the holonomy group |T|. Each single image 
sees mult(#) neighbours lying at distance 9 from it, so that this distance appears 
|r| mult(#)/2 times (the division by two compensates for the fact that we counted the 
distance from image A to image B as well as the distance from image B to image A). 
We conclude that if we have N sources in the fundamental polyhedron, the amplitude 
of the spike located at x = is 



As an example, let us consider the case of the cyclic group Z e consisting of six 
Clifford translations through distances (— 27r/3, — 7r/3, 0, 7r/3, 27r/3, 7r). These Clifford 
translations yield spikes of multiplicity 1 at x = and x = tt, and spikes of multiplicity 
2 at x = tt/3 and \ = 2n/3 7 as shown in table 2. Applying the above formula to the 
spike at x — tt/3, in the case of N = 300 distinct sources in the fundamental domain, 
we expect the spike to reach a height \ N \T\ mult(0) = ±(100)(6)(2) = 1800 above 
the background distribution, in agreement with computer simulations (see figure 13). 



T > {Xa,Xl) 




0(271" - 2xa - Xi) x [ sin2xa + 7T - Xa - Xi/2 
- cosxa sec(x;/2) sin(x a - Xj/2)]}, 



(22) 
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1/2 1 


multiplicity 


1 6 1 



Table 1. Position and the multiplicity of the spikes for the binary dihedral group 
D|. Compare with figure 14. 





1/3 2/3 1 


multiplicity 


12 2 1 



Table 2. Position and the multiplicity of the spikes for the binary cyclic group 
Z(j. Compare with figure 14. 



xA 


1/3 1/2 2/3 1 


multiplicity 


1 2 6 2 1 



Table 3. Position and the multiplicity of the spikes for X Zz- As expected, 
we find that it can be obtained from the data for the binary dihedral group 
(table 1) and Z3 X {±1} = Z§ (table 2). Compare with figure 14. 



xA 


1/5 1/3 2/5 1/2 3/5 2/3 4/5 1 


multiplicity 


1 12 20 12 30 12 20 12 1 



Table 4. Position and the multiplicity of the spikes for the binary icosahcdral 
group I*. Compare with figure 15. 

In tables 1 to 4, we give the translation distances |0| and the multiplicities mult(0) 
for the binary dihedral group the cyclic group Z§, the product D 1 ^ x Z3, and the 
binary icosahedral group /*, respectively. They correspond to the PSH spectra shown 
in figures 14 and 16. Note that the spectrum for D\ x Z 3 (table 3) is obtained from 
those of the binary dihedral group (table 1) and Z 3 x {±1} = Z & (table 2). 

5.4- PSH numerical Simulations 

The previous section presented the theoretical expectations for PSH in multiply 
connected spherical spaces. Here we carry out numerical simulations confirming those 
expectations. Later in this section we will also take into account the approximate 
flatness of the observable universe, which implies that we are seeing only a small 
part of the covering space, and therefore can observe spikes only at small comoving 
distances \- 

In figure 13 we draw the histogram for the lens space L(6, 1) and we check that 
our geometrical expectations about the positions and the heights of the spikes are 
satisfied. 

In figures 14 to 16 we present various simulations. We have chosen to draw 
the normalised PSH in the covering space, i.e. the PSH divided by the number of 
pairs and the width of the bin, and to show in grey regions the part of the PSH 
that is observable if the redshift cut-off is respectively z c = 1,3, 1000, which roughly 
corresponds to catalogs of galaxies, quasars and the CMB. 

We start by showing in figure 14 the "additivity" of the spectrum by considering 
the group D% x Z3 for which the positions of the spikes in its PSH can be found 
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from the PSH of the binary dihedral group D'^ and of Z 3 x {±1} = Z 6 . Indeed, the 
amplitude of the spikes is smaller since the same number of pairs has to contribute to 
more spikes. 

We then show in figures 15 and 16 different groups of the same order, chosen to be 
120. The aim here is first to show that the number of spikes is not directly related to 
the order of the group and that the order of the cyclic factor, if any, is more important. 
While increasing this order, the number of spikes grows but their amplitudes diminish. 
To get a physical understanding of this effect, we show in figure 17 the view for an 
observer inside the manifold for the binary octahedral group O*, for the cyclic group 
Zn, and then for their product. 

Finally, we applied our calculations by applying the PSH to real data, namely 
a catalog of about 900 Abell and ACO clusters with published rcdshifts. The depth 
of the catalog is z cut — 0.26, corresponding to 730/i _1 Mpc in a spherical universe 
(see [14] for a more detailed description of this catalog). The PSH exhibits no spike; 
this gives a constraint on the minimal size of spherical space, which corresponds to a 
maximum order of the cyclic factor of the holonomy group of about 200. 
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Figure 13. The PSH for the lens space L(6, 1) generated by the cyclic group 
Zq with N = 300 sources in the fundamental domain. The height of the peaks 
at x = n /3 and x = 27r/3 is of order 4200 while the background distribution is 
of order 2400, giving the peaks a height above the background of about 1800, in 
agreement with our theoretical analysis. 

5.5. Robustness of PSH 

In [20] we addressed the question of the stability of crystallograpic methods, i.e. to 
what extent the topological signal is robust for less than perfect data. We listed the 
various sources of observational uncertainties in catalogs of cosmic objects as: 

(A) the errors in the positions of observed objects, namely 

(Al) the uncertainty in the determination of the redshifts due to spectroscopic 
imprecision, 

(A2) the uncertainty in the position due to peculiar velocities of objects, 
(A3) the uncertainty in the cosmological parameters, which induces an error in 
the determination of the radial distance, 
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Distance/ft 

Figure 14. The PSH spectra for (upper left) the binary dihedral group DJ, 
(upper right) the cyclic group Zg and (bottom) the group X Z3. One can 
trace the spikes from each subgroup and the black line represents the analytic 
distribution in a simply-connected universe. 

(A4) the angular displacement due to gravitational lensing by large scale structure, 
(B) the incompleteness of the catalog, namely 
(Bl) selection effects, 

(B2) the partial coverage of the celestial sphere. 

We showed that the PSH method was robust to observational imprecisions, but 
able to detect only topologies whose holonomy groups contain Clifford translations, 
and we performed numerical calculations in Euclidean topologies. Fortunately the 
shortest translations in spherical universes are typically Clifford translations (even 
though more distant translations might not be), so the PSH is well suited to detecting 
spherical topologies. In the present work we check the robustness of PSH in spherical 
topologies. 

Let us consider various spherical spaces of the same order 120 - namely the 
lens space L(120, 1), the binary dihedral space D^ and the Poincare space I* - and 
calculate the corresponding PSH's. The density parameters are fixed to f2 TOo = 0.35 
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Distance/^ Distance/71 

Figure 15. The PSH spectra for two spherical spaces of order 120 but with 
no cyclic factor: (left) the binary dihedral group D| and (right) the binary 
icosahedral group I* 



Table 5. Values of the limit percentage pi of rejection above which the PSH 
spikes disappear, as a function of the redshift cut-off for the Poincare space I*. 



^cut 


1 


3 


> 5 


Pi 


No signal 


80% 


90% 



and il\ = 0.75. In the runs, the number of objects in the catalog is kept constant. 
We examine separately the effects of errors in position due to redshift uncertainty Az, 
and the effects of catalog incompleteness. Each of these effects will contribute to spoil 
the sharpness of the topological signal. For a given depth of the catalog, namely a 
redshift cut-off z cut , we perform the runs to look for the critical value of the error at 
which the topological signal fades out. 

Figure 18 gives the critical redshift error Azi above which the topological spikes 
disappear for the spherical space Z120 with Q — 0.35, Cl\o = 0.75. 

Next, we simulate an incomplete catalog where we randomly throw out p% of the 
objects from the ideal catalog. For the cyclic group .Z120 and the binary dihedral group 
D| , the topological signal is destroyed only for a very large rejection percentage (above 
90%). For the less favorable case of the Poincare group I*, the results are summarized 
in table 5. Thus our calculations confirm that the PSH method is perfectly robust for 
all spherical topologies containing Clifford translations. 

As we have seen the PSH method applies for most of the spherical manifolds. 
Nevertheless when the order of the group or of one of its cyclic subgroups is too 
high then the spikes are numerous and have a small amplitude. Thus they may be 
very difficult to detect individually. This the case for instance for some L(j>, q) as 
well as for linked action manifolds. In that case the CCP method, which gathers the 
topological signal into a single index, is more suitable. Again the topological signal is 
rather insensitive to reasonable observational errors as soon as the underlying geometry 
contains enough Clifford translations. 
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Figure 16. The PSH spectra for three spaces of order 120 with different cyclic 
factors: (upper left) Dg X (upper right) Dg X Z^ 2 and (bottom) Zi2o- As 
explained in the text, the order of the group being fixed, the higher the order of 
the cyclic subgroup the larger the number of spikes. However the amplitude of 
these spikes is smaller. 



6. Conclusion and Perspectives 

In this article, we have investigated the possible topologies of a locally spherical 
universe in the framework of Friedmann-Lemaitre spacetimes. We have given the 
first primer of the classification of three-dimensional spherical spaceforms, including 
the constructions of these spaces. 

We have determined the topologies which are likely to be detectable in three- 
dimensional catalogs of cosmic objects using crystallographic methods, as a function 
of the cosmological paramaters and the depth of the survey. The expected form of the 
Pair Separation Histogram is predicted, including both the background distribution 
and the location and height of the spikes. We have performed computer simulations of 
PSH in various spherical spaces to check our predictions. The stability of the method 
with respect to observational uncertainties in real data was also proved. 

Such a complete and exhaustive investigation of the geometrical properties of 
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Figure 17. To have an intuitive understanding of which subgroups give rise to 
spikes and how, it is fruitful to have a look from the inside. In the upper left panel, 
we see the view inside the quotient of the 3— sphere S 3 by the binary octahedral 
group O*. The fundamental domain is a truncated cube, 48 copies of which tile 
S 3 . The tiling's 1-dimensional edges are shown in the figure. In the upper left 
panel, we see the view in 1/(17, 1). All 17 translates of the Earth align along a 
Clifford parallel. When mixing both groups to get O* X Z17 in the bottom panel, 
we see the simultaneous effects of the cyclic factor Z17, which generates lines of 
images, and the binary octahedral factor, which translates one line of images to 
another. This illustrates how important the order of the cyclic group is, because 
it alone determines the distance to our nearest translate. 

Note: More distant images of the Earth are always dimmer than closer images 
(because of artificial "fog" ) , but the apparent size of an image decreases only until 
the images reaches a distance of it/ 2, after which more distant images appear 
larger because the light from them follows geodesies that reconverge in S 3 . As an 
image approaches the antipodal point, at a distance of 7T, it fills the whole sky. 



spherical spaces will be useful for further studies dealing with spherical topologies, 
including two-dimensional methods using CMB data. We plan to investigate the 
applicability of the circle-matching method [40] for spherical topologies. 




Figure 18. Plot of Az; as a function of the depth of the catalog, for the spherical 
space Z120 with Qq = 0.35, Qao = 0.75, using the PSH method. 



Appendix A. Quaternions 

The quaternions are a four-dimensional generalization of the familiar complex 
numbers. While the complex numbers have a single imaginary quantity i satisfying 
i 2 = — 1, the quaternions have three imaginary quantities i, j, and k satisfying 

i 2 = j 2 = k 2 = -1 (A.l) 
which anti-commute 

{i,j} = 0, {j,k} = 0, {k,i} = (A.2) 
and are subject to the multiplication rules 

[i,j] = 2k, [j,k] = 2i, [k,i] = 2j, (A.3) 

and 

[i,l] = 0, [j, 1] = 0, [k,l]=0, (A.4) 

where [] and {} are the usual commutation and anti-commutation symbols. 
Geometrically, the set of all quaternions 

q = al + M + cj + dk, (a, b, c, d) e R* (A. 5) 

defines four-dimensional Euclidean space, and the set of all unit length quaternions, 
that is, all quaternions al + bi + cj + dk satisfying a 2 + b 2 + c 2 + d 2 = 1, defines the 
3-sphere. The quaternions are associative, (qr)s = q(rs), even though they are not 
commutative. Introducing the conjugate quaternion q of q by 

q = al — bi — cj — dk (A. 6) 

and the modulus of q by 



|q| = Vqq= \/a 2 + b 2 + c 2 + d 2 (A.7) 
a unit quaternion satisfies 

|q| = l^q = q- 1 . (A.8) 
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The identity quaternion 1 is fundamentally different from the purely imaginary 
quaternions i, j and k, but among the unit length purely imaginary quaternions 
bi + cj + dk there is nothing special about the basis quaternions i, j and k. Any 
other orthonormal basis of purely imaginary quaternions would serve equally well. 
Lemma Al: Quaternion change of basis 

Let M be a 3 x 3 orthogonal matrix, and define 
i' = M n i + M 12 j + M 13 k 
j' = M 2 ii + M 22 j + M 23 k 

k' = M 3 ii + M 32 j + M 33 k (A.9) 
then i', j' and k' satisfy the quaternion relations 
i' 2 = j' 2 = k' 2 = 1 
i'j' = k' = -j'i' 
j'k' = i' = -k'j' 

k'i' =j' = -i'k' (A.10) 

Lemma Al says that an arbitrary purely imaginary quaternion oi+cj+dk may, by 
change of basis, be written as b'i' . If the purely imaginary quaternion M + cj + dk has 
unit length, it may be written even more simply as i'. A not-necessarily-imaginary 
quaternion al + bi + cj + dk may be transformed to a'l + b'i'. If it has unit length it 
may be written as cos0 1 + sin0 i' for some 0. 

The unit length quaternions, which we continue to visualize as the 3-sphere, may 
act on themselves by conjugation or by left or right multiplication. 
Proposition A2: Conjugation by quaternions 

Let q be a unit length quaternion. According to the preceding discussion, we may 
choose a basis {l,i',j',k'} such that q = cos0 1 + sin0 i' for some 0. It is easy to 
compute how q acts by conjugation on the basis {1, i', j', k'}: 

(cos 1 + sin i') 1 (cos 01 -sin i') = 1 (A.ll) 
(cos 9 1 + sin 9 i')i' (cos 91 -sin 9 i') = i' (A.12) 
(cos 6*1 + sin0i')j' (cos 01 - sin0i') = cos 20 j' + sin 20 k' (A.13) 
(cos 9 1 + sin 9 i') k' (cos 9 1 - sin0 i') = - sin 29 j' + cos 29 k' (A.14) 
Conjugation by any quaternion fixes 1 ("the north pole"), so the action is always 
confined to the "equatorial 2-sphere" spanned by {i',j',k'}. Within the equatorial 
2-sphere, conjugation by the particular quaternion cos 91 + sin 9 i' rotates about the 
i' axis through an angle 29. Unlike the preceding action by conjugation, which always 
has fixed points, when the quaternions act by left or right multiplication they never 
have fixed points. 

Proposition A3: Left multiplication by quaternions 

Let q be a unit length quaternion. Choose a basis {l,i',j',k'} such that q = 
cos0 1 + sin0 i' for some 9. It is easy to compute how q acts by left multiplication on 



the basis {1, i', j', k'}. 

(cos 91 + sin0i')l = cos 9 1 + sin 9 i' (A.15) 

(cos 01 + sin0i')i' = - sin0 1 + cos0i' (A.16) 

(cos 01 + sin0i')j' = cos0 j' + sin k' (A.17) 

(cos 1 + sin i')k' = - sin j' + cos k'. (A.18) 
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We see that left multiplication rotates 1 towards i' while simultaneously rotating j' 
towards k'. The result is a screw motion along so-called Clifford parallels. The Clifford 
parallels are geodesies, and they are homogeneous in the sense that there is a parallel- 
preserving isometry of 5 3 taking any one of them to any other (see section 3.1). Action 
by right multiplication is similar, but yields left-handed Clifford parallels instead of 
right-handed ones. 



Appendix B. Matrices 



This appendix defines each finite subgroup of 5*0(3) as an explicit set of rotations. 

Lifting a subgroup of SO (3) to the corresponding binary subgroup of S 3, is easy. 
Proposition A. 2, together with the change of basis principle in Proposition A.l, implies 
that a rotation through an angle 9 about an axis (x, y, z) is realized by both the unit 
length quaternion 

6 1 9 

q = cos - 1 H ■j===== sin - (x i + y j + z k) (B.l) 



y 



= sin ■ 

2 + z 2 2 



and by its negative. This one-to-two mapping yields a binary subgroup of S 3 that is 
exactly twice as big as the original subgroup of 5*0(3) . 

Converting from a quaternion in S 3 to a matrix in 50(4) is also easy. When a 
quaternion q = al + bi + cj + dk acts on the left its effect on the basis vectors 1, i, j, 
and k is, respectively, 

(al + bi + cj + dk) 1 = al + bi + cj + dk 
(al + bi + cj + dk) i = —61 + ai + dj — ck 
(al + bi + cj + dk) j = -cl - di + aj + 6k 

(al + bi + cj + dk) k = -dl + ci - 6j + ak (B.2) 
Thus, relative to the orthonormal basis {1, i, j, k}, the matrix for this action is 



/ 



Mieft(q) 



c 



-b 
a 
d 

— c 



— c 
-d 
a 
6 



-d \ 

c 

-6 

a 



(B.3) 



Similarly, when the same quaternion q = al 
resulting isometry has matrix 



M righ t(q) = 



/ 


a 


-b 


— c 


-d\ 




6 


a 


d 


— c 




c 


-d 


a 


6 


V 


d 


c 


-b 


a J 



bi + cj + dk acts on the right, the 



(B.4) 



In the case of double action and linked action manifolds (§4.2 and §4.3), one 
quaternion r acts on the left (giving a right-handed Clifford translation) while a 
different quaternion / acts on the right (giving a left-handed Clifford translation). 
Their combined action is given by the product of their respective matrices. The order 
in which we multiply the two matrices doesn't matter because right and left handed 
Clifford translations always commute. 

The following realizations of the finite subgroups of 50(3) are unique up to an 
orthonormal change of coordinates. 
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Appendix B.l. Cyclic Groups Z n 

Each cyclic group Z n consists of rotations about the axis (0,0,1) through an angle 
27rfc/n, for < k < n. This defines n distinct rotations. 

Appendix B.2. Dihedral Groups D m 

The dihedral group D m contains the m rotations of the cyclic group Z m about the 
axis (0,0,1) as well as m half turns about the axes (cos ^ , sin ^ , 0) for < k < m. 
This defines 2m distinct rotations. 

Appendix B.3. The Tetrahedral Group T 
The tetrahedral group (Figure Bl) consists of 

• the identity, 

• order 2 rotations about the midpoints of the tetrahedron's edges, realized as it 
rotations about the three axes (1,0,0), (0,1,0), and (0,0,1), and 

• order 3 rotations about the tetrahedron's vertices and face centers, realized as 
counterclockwise 27r/3 rotations about the eight axes (±1,±1,±1). 

This defines a total of 1 + 3 + 8 = 12 = |T| distinct rotations. Applying equation (B.l) 
yields the 24 quaternions of the binary tetrahedral group T* . 




Figure Bl. Visualisation of the symmetry group T of the tetrahedron. 

Appendix B.4- The Octahedral Group O 
The octahedral group consists of 

• the identity, 

• order 2 rotations about the midpoints of the octahedron's edges, realized as 7r 
rotations about the six axes (±1,1,0), (0, ±1,1), and (1,0, ±1), and 

• order 2 rotations about the octahedron's vertices, realized as 7r rotations about 
the three axes (1,0,0), (0,1,0), and (0,0,1), and 

• order 3 rotations about the centers of the octahedron's faces, realized as 
counterclockwise 27r/3 rotations about the eight axes (±1,±1,±1). 
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• order 4 rotations about the octahedron's vertices, realized as counterclockwise 
7r/2 rotations about the six axes (±1,0,0), (0, ±1,0), and (0,0, ±1). 

This defines a total of 1 + 6 + 3 + 8 + 6 = 24 = \0\ distinct rotations. Applying 
equation (B.l) yields the 48 quaternions of the binary octahedral group O* . Note 
that the octahedral group contains the tetrahedral group within it. 

Appendix B.5. Icosahedral Group I 

The icosahedron's twelve vertices being at (±0, ±1,0), (O,±0, ±1), and (±1,0, ±0), 
where = (V5 — l)/2 is the golden ratio, the icosahedral group consists of 

• the identity, 

• order 2 rotations about the midpoints of the icosahedron's edges, realized as 
7r rotations about the fifteen axes (1,0,0), (0,1,0), (0,0,1), (±0,±(0 + 1),1), 
(1, ±0, ±(0 + 1)), and (±(0 + 1), 1, ±0), and 

• order 3 rotations about the centers of the icosahedron's faces, realized as 
counterclockwise 27r/3 rotations about the twenty axes (±1,±1,±1), (O,±(0 + 
2),±1), (±1,0, ±(0 + 2)), (±(0 + 2), ±1,0), and 

• order 5 rotations about the icosahedron's vertices, realized as both counterclock- 
wise 2tt/5 rotations and counterclockwise 47r/5 rotations about the twelve axes 
(±0,±1,O), (O,±0,±l), and (±1,0, ±0). 

This defines a total of 1 + 15 + 20 + 24 = 60 = |/| distinct rotations. Applying 
equation (B.l) yields the 120 quaternions of the binary icosahedral group I*. Note 
that the icosahedral group also contains the tetrahedral group within it. 
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